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Abstract 

It is the intention of this paper to rigorously clarify the role of the occupation 
numbers in the current practical applications of the density functional formal- 
ism. In these calculations one has to decide how to distribute a given, fixed 
number of electrons over a set of single-particle orbitals. The conventional 
choice is to have single-particle orbitals below the Fermi level completely occu- 
pied (with possible fractional occupations at the Fermi level) and the orbitals 
above the Fermi level empty. Although there is a certain confusion in the lit- 
erature why this choice is superior to any others, the general belief is that it 
can be justified by treating the occupation numbers as variational parameters 
and then applying Janak’s theorem or similar reasoning. We demonstrate 
that there is a serious flaw in those arguments, mainly the kinetic energy 
functional and therefore the exchange-correlation functional are not differ- 
entiable with respect to the density for arbitrary occupation numbers. It is 
rigorously shown that in the present context of the density functional calcu- 
lations there is no freedom to vary the occupation numbers. The occupation 

numbers cannot be considered as variational parameters. 
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I. INTRODUCTION. 



Predicting properties of real materials from first principles has been an important goal of 
solid state physics from the early days, going back to Thomas, Fermi and others. During the 
past few decades, density functional theoryEhi (DFT) within the local density approximation 
(LDA) has been quite successful in this regard. This method has no adjustable parameters 
and the only significant approximation is the local density approximation. However this 
method also relies heavily on the variational nature of the ground state total energy and 
this paper is concerned with concepts and misconcepts related to this variational property. 

Nonvariational nature of the energy functional with variable occupation numbers has 
been noticed recently by Weinert and Davenporti. However we will rigorously prove that 
the differentiability of the energy functional with respect to the density depends critically 
on the choice of the occupation numbers and that there is no freedom to vary these without 
sacrificing the differentiability of the energy functional, i.e. a variational principle does not 
exist for any arbitrary choice of occupation numbers within the current framework of the 
density functional calculations. As we show later this means that Janak’s theorem cannot 
exist within the current framework of density functional calculations. To our knowledge the 
variational nature of the occupation numbers has never been studied in this fashion. The 
significance and relevance of the above statements should be judged in the context of the 
references made to this problem over the past two decadesIHl. 

Original Hohenberg-Kohn-Sham theoryB! (HKS) rests on the assumption that the ground 
state density of the interacting system is simultaneously a ground state density of the non- 
interacting system with suitably chosen external potential. The ground state for the non- 
interacting system can always be obtained by completely filling the lowest single-particle 
orbitals. The occupation numbers in this approach are well defined (except for the case 
when the topmost level is degenerate) and any other choice would be in direct violation 
with the original HKS assumption. In order to make any rigorous statements about this 
particular choice for the occupation numbers one has to define more general energy func- 
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tional with the domain that includes densities which are not the ground state densities of 
any noninteracting system. Such energy functional can be defined by using the constrained 
search approach proposed by Levyi . 



II. CONSTRAINED SEARCH APPROACH. 

Constrained search approach! defines the energy functional in a very simple and a phys- 
ically appealing way. The ground state of any many-electron system can be found as 

E 0 = mm(ty\H\ty). (1) 

Here the minimization is performed over all allowable antisymmetric wave functions. To 
include mixed states as well as pure states one can extend the minimization from wave 
functions to all allowable density matrices: 

Eq = ruin Tr{YH}. (2) 

f 

where 

r = 'Edi\v i )(Vi\, (3) 

di = d* > 0, = 1- 

The minimization can be divided into two steps: 

£0 = min{minTr{M}}. (4) 

n f — m 

The first minimization is performed over all allowable density matrices that lead to some 
fixed density n(r). The result of this minimization is the functional of the density only: 

El[u ] = mmTr{f H}. (5) 

f— >n 

According to eq. (4) once we have found this functional, the ground state energy as well as 
the ground state density can be found by minimizing eq. (5) with respect to the density, i.e. 
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Eq = min El \n] . 

n 



( 6 ) 



The Hamiltonian for the many-electron system in general can be written as 



H = T + W + V, 



(7) 



where 



is the kinetic energy term, 



T = — - J 'il>\r)'V 2, ip(r)dr 



W — — J -'0'^(r)'i/y*'(r / )'0(r , )'i/t(r)(ir 



is the electron-electron interaction term, 



( 8 ) 



(9) 



V — n(r)'0^(r)'0(r)hr 



( 10 ) 



is the external field interaction term. According to eq.(5) this leads to 



E L [n] = F L [n] + J v(r)n(r)dr. 



( 11 ) 



All the complexities of many-electron system are now hidden in the functional Fl [n] : 



F L [n] = min TrjfT + fw}. 



r^n 



( 12 ) 



It was shown by Lieb0 that Fl[u\ and therefore E Iy [n] are convex functionals. This means 
that Ei[n ] has no extrema above its absolute minimum0. 

The following important theorem by H.Englisch and R.Englisch@ is central to 
elusion and deals with the issue of the differentiability of Fl [n ] . 

Theorem 1: The functional 



our con- 



Fl [n] = minTr{fT + f IE} 

f 1 — >n 



is differentiable nowhere else but on the following set of densities: 
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n(r) = ^ A fc n fc (r), 

k 

with J2k — 1, Afe > 0, and nk{ r) = / |^(r, r 2 , r/v)| 2 dr 2 ...<irjv ; where ifk are the or- 
thonormal, degenerate ground states of the Hamiltonian given by (7). Moreover the choice 
for the external potential is unique (up to some additive constant) for a given density. 

As a corollary to the above theorem one can prove the following statement: 
Proposition 1: The functional 



Tj[n] 



mm 

Yji 



y 



fi / (r) 




is differentiable nowhere else but on the following set of densities: 



(13) 



where 



n(r ) = ^/i|<Mr)| 2 , 



— 7T + V ( r ) <M r) = r), 



and Yffi = N, 

1 for 6j < e F 

fi — * 0 for 6i > e F (14) 

x for €i = e F , where 0 < x < 1 

The proof of the above statement follows from the fact that ensemble-N-representable first 
order density matrix (i.e derivable from N-particle density matrix T) can always be written 



7 <r,r') = y/^*(r')*(r), 



(15) 



where 0 < f < 1 , fi — AT and <f>f s are orthonormal. Since the expectation value of the 
kinetic energy depends only on the first order density matrix it follows that@ 



Tj[n ] = min Tr {TT}. 

T— >n 
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Once this equivalence has been established, we can now apply Theorem 1 (with W set 
to zero) to the functional Tj[n]. In this case if>k s are the ground states of the following 
Hamiltonian 

N i N 

i= 1 Z i = 1 

and can be represented as single determinants built from N lowest single particle orbitals 

</>(r): 

(-^V 2 + v(r))0j(r) = 6i<t)i{ r) ei < e 2 < e 3 < .... 

The only case where we have degenerate many-particle ground state is when there is a 
degenerate single particle level (k < N ) with degeneracyS m > N — k + 1. This energy 
level 6k is the Fermi level (e& = e/). Then 

M r ) = /|'0fe(r,r 2 ,...,r A r)| 2 dr 2 ...dr J v = ^ \M r )\ 2 + Y l^( r )| 2 - 

e ij c =e F 

The densities where the functional Tj[n] is differentiable are given by: 



where 



™( r ) = H^fc( r ) = E A fc H l^( r )| 2 + 5Z A fc H l<M r )| 2 

/c k €i<6F k e i k =e F 

= Y l<M r )l 2 + fi\M r ) I 2 , 

£;<£F £i=£F 



k 



(16) 



III. JANAK’S THEOREM. 



It is obvious that the energy functional (11) is of little help when it comes to practical 
applications. Further simplifications are needed. Following the HKS strategy one usually 
IiartitionJzBO the functional Fj j [ n\ in the following fashion: 

1 



F L \n\ = Tj[n\ + 



n(rWr') , , , „ r , 

-1-drdr + E Lxc [n\. 



r — r 



(17) 
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This results in the following energy functional 
E[n } = Tj[n } + J v(r)n(r)dr + ; 



n(r)n(r') , , , ^ r , 

-1-drdr + E Lxc [n\ 



r — r 



( 18 ) 



If the above functional was differentiable for any choice of the occupation numbers, its 
subsequent minimization with respect to density would lead@’0E! to the following results: 



-~2 + v *ff( T ) J <M r ) = et0*(r), 



(19) 



where 



v eff(r) = v(r) + J 



n ( r ') dr , + SE Lxc [n] 



r — r 



5n( r) 



n ( r ) = 



and 



SE[n] 



= e*. 



( 20 ) 



The last result was first obtained by Janak@ before the advent of the constrained search 
formulation and is known as Janak’s theorem. This theorem is usually regarded as a justi- 
fication for occupying only the single particle orbitals below the Fermi level-!. 

Nevertheless, the flaw in the above derivation is obvious. According to Proposition 1, 
once the partition (17) has been performed we have immediate restrictions on the occupa- 
tion numbers. Any choice other than the conventional one, makes functionals Tj[n ] and 
Elxc\p\ (since El xc [vi\ implicitly contains the kinetic energy term Tj[n ]) nonclifferentiable 
with respect to density. Therefore we cannot vary the energy functional (18) for some arbi- 
trary choice of the of the occupation numbers {/*}. Hence equation (19) is not valid for the 
occupation numbers other than the ones specified in Proposition 1. Obviously, in this case 
eq. (20) is not valid since we no longer have the freedom to choose arbitrary occupation 
numbers. In particular this means that there is no room for the Janak’s theorem within the 
current framework of density functional calculations. 

The question now arises how to occupy the topmost degenerate levels (if any). As 
previously mentioned, the functional El[ti] is convex and there are no extrema above its 
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absolute minimum. Therefore any fractional occupations at the Fermi level that conserve 
the degeneracy (i.e. the functionals Tj[n ] and El xc [ii\ remain differentiable) and leave eq. 
( 17 ) intact (i.e. deliver the extremum of the functional Ei [n] ) should give the same total 
energy0. This of course may not be true in practice because of our approximations to the 
exchange-correlation energy. 



IV. CONCLUSION. 

The functionals Tj[n\ and E Lxc [n ] appear as a direct consequence of the attempt to 
map the interacting system into the noninteracting one. We have demonstrated that these 
functionals are differentiable only on a certain domain of densities. It is only there that 
we are allowed to vary the energy functional. It is shown that in order for density to 
belong to this domain there is a certain restriction on the occupation numbers, i.e. single- 
particle orbitals below the Fermi level are completely occupied (with possible fractional 
occupations at the Fermi level) and the orbitals above the Fermi level are empty. Any other 
choice of the occupation numbers will make the functionals Tj[n\ and El xc [ti\ lose their 
differentiability with respect to density - a property which is imperative for the variational 
principle. Note that even though the approximate forms for the exchange-correlation energy 
currently being used might appear differentiable with respect to density regardless of the 
choice of the occupation numbers, one can hardly claim that these approximations will hold 
for an arbitrary set of the occupation numbers. 

The main conclusion is that the occupation numbers cannnot be considered as variational 
parameters. It means that there is no apriory justification for using the conventional set of 
the occupation numbers in the current density functional calculations. The only reason we 
resort to this set of the occupation numbers is because any other choice would make the 
variational principle inapplicable. 
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